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Lattice QCD in rotating frames 
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We formulate lattice QCD in rotating frames to study the physics of QCD matter under rotation. 
We construct the lattice QCD action with the rotational metric and apply it to the Monte Carlo 
simulation. As the first application, we calculate the angular momenta of gluons and quarks in the 
rotating QCD vacuum. This new framework is useful to analyze various rotation-related phenomena 
in QCD. 
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Introduction. — Quantum Chromodynamics (QCD), 
the fundamental theory of the strong interaction, has 
been resulted in incredibly diverse phenomena. Among 
them, QCD matter under rotation is of particular in- 
terest. For example, the cores of rapidly rotating com- 
pact stars arc expected to be in the superfluid or su- 
perconducting phase and arc threaded with topological 
vortices [l[, which may play an important role in pul- 
sar glitches m. In ultrarelativistic heavy- ion collisions 
performed at Relativistic Heavy Ion Collider (RHIC) or 
Large Hadron Collider (LHC), the created quark-gluon 
plasma should have a finite angular momentum, espe- 
cially in off-central collisions. This would result in inter- 
esting phenomena such as the chiral vortical effect Q or 
the Kelvin- Helmholtz instability Also in low-energy 
nuclear physics, rotation generates characteristic states 
of nucleus, e.g., rotational modes and high-spin states 
0. Thus, the analysis of the QCD matter under rotating 
environment is an important phenomenological problem. 

However, first principle calculations for the rotating 
QCD matter has been elusive. In numerical simulations 
under equilibrium conditions, such as the lattice QCD 
simulation, it is not straightforward to generate a rotat- 
ing state. The rotating state cannot be realized by adding 
external fields, since the circulating velocity field, which 
characterizes the rotation, is the matter field itself. One 
possible way is to rotate the reference frame. A rest state 
in a rotating frame is equivalent to a rotating state in a 
rest frame. By performing the simulations in rotating 
frames, we can study the properties of the matter under 
rotation. 

In this Letter, we formulate lattice QCD in rotating 
frames. This is one practical application of lattice QCD 
in a curved space-time, i.e., in general relativity. Lattice 
QCD in moving frames has been discussed only in the 
case of the Lorentz boost, i.e., in special relativity 
Lattice gauge theory in a curved space-time has been 
discussed in the context of lattice quantum gravity • 

In the continuum. — In the standard lattice QCD sim- 
ulation, we numerically perform the path integral with 
the Euclidean metric gij = diag(l, 1, 1, 1). To formu- 



late rotation in the Euclidean space-time, we need two 
operations: the Wick rotation r = —it and the spatial 
rotation 9 = di-cst — ^t- There are two possibilities for 
the order of these two operations: the "Minkowskian" 
rotation ft = —dO/dt, which is defined as the spatial 
rotation before the Wick rotation, and the "Euclidean" 
rotation il ~ —dO/dr, which is defined as the spatial 
rotation after the Wick rotation. As is shown below, 
there is the sign problem in the Minkowskian rotation 
and no sign problem in the Euclidean rotation. (This re- 
lation is similar to the relation between the Minkowskian 
and Euclidean electric fields [sj-) As long as the analytic 
continuation to the original Minkowski space-time is val- 
idated, these two rotations produce the same end result. 
In this Letter, we apply the Euclidean rotation to the 
lattice QCD simulation. In the following equations, the 
angular velocity of the Euclidean rotation is denoted by 
51. The corresponding equations with the Minkowskian 
rotation is obtained by replacing 51 — !■ i51. 

Wc choose the Cartesian coordinate a;^ = {x,y,z,T). 
The metric tensor of a rotating frame is 
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Here r = ^ x"^ Ar -ip- is the distance from the rotation 
axis. In the rotating frame, an observer feels the local 
Lorentz boost with the Lorentz factor + r^il^. In 

the Minkowskian rotation, for a well-defined coordinate 
patch, we must impose the condition rJl < 1 = c (the 
light velocity), which means that the local velocity must 
be smaller than the light velocity. In the Euclidean rota- 
tion, there is no such coordinate singularity but the local 
velocity should be small for the analytic continuation. 

We start with the gluon and quark actions in a general 
curved space-time, 

Sg = f d^x^dctgo^p -^^g^^gP^ivF^pF,,, (2) 
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d^x^Jdeig^p i>[l^{D^ - F^) + m]?/>. (3) 
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The covariant Dirac operator is constructed from the 
SU(-/Vc) covariant derivative = — iA^ and the 
spinor affine connection F^. The connection is defined 
as 



The lattice is rotated around the z axis. In the x and y 
directions, we take the Dirichlet boundary conditions. In 
the z and r directions, we take boundary conditions in 
the same manner as the usual lattice simulation. 



7^7'), 



(4) 

(5) 
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where ef is the vierbein and F^^ is the Christoffel sym- 
bol. The Greek and Latin indices refer to the coordinate 
and tangent spaces, respectively. 

By substituting the rotational metric ([1]), the gluon 
action is 

1 



Sg = / d^x —tT[{l + r'^')F,yF,y 
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(7) 



(1 + y^n^)F,,F,, + (1 + x^n^)Fy,Fy, 

^y^FxyFy-j- '2x^Fyx FxT 

2ynF^,F,r - 2xnFy,F,r + 2xyn^F.,,F,y]. 



As an effect of rotation, the gluon action includes the 
0(0) terms which break parity and time- reversal sym- 
metry and the O(ri^) terms which do not break them. 

The covariant Dirac operator depends on the choice of 
the vierbein. We choose the vierbein 



^X — _y — ^Z — _T — 



1, 



64 = -yQ, el xn, 



(8) 
(9) 



and 



for other components. In this choice, the 



quark action is 



d x ip 



^-D^ + ^yDy + YD, 

(10) 



Y Dr + in- 



The gamma matrices in the rotating frame are given as 



, i.e. 
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7 -y^i , 
-f^ + xn-f\ 

7^ 



(11) 

(12) 
(13) 
(14) 



As a result of rotation, the Dirac operator includes the 
orbit- rotation coupling term ^'^U,{xDy — yD^) and the 
spin-rotation coupling term i^'^fla^'^ /2. 

On the lattice. — We discretize the continuum actions 
([7]) and (|10|) on the hypercubic lattice. The schematic 
figure is shown in Fig.[T] The lattice spacing is a and the 
total number of the lattice sites is x Ny x x Nr- 




FIG. 1: Rotating lattice. 



On the lattice, the gluon field strength is constructed 
from the gauge invariant loops of the link variables Uf^{x). 
The squared terms, e.g., F^yF^y, are constructed from 
the" plaquette" and the non-squared terms, e^., F^yFy^, 
arc constructed from the "chair- type" loop [9|. The lat- 
tice gluon action is 
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+ {l + y^n')[l 



-Retr[/^ 



+ {l+x^n') [ l-—RetTUy, 



3 - —Re tr ([/^^ + Uyr + U,r) 



(15) 
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■Retr(yr2I/j;j,T- — xftVy, 



y^VoozT - xflVyzr + XyQ'^Vooz 



The bare lattice coupling is /3 = 2Nc/ cjym.- ^'^^ ^ local 
definition of the lattice field strength, we take the clover- 
type average of four plaquettes as 
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(16) 



In a flat space-time, this average is redundant because 
it can be absorbed into the summation in the ac- 
tion. However, in a curved space-time, this average is 
important because the coefficients depend on the space- 
time coordinate. Similarly, we take the (anti-)symmetric 
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average of eight chair-type loops as 
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(17) 



For the lattice quark action, we adopt the Wilson fermion 
with the gamma matrices in the rotating frame. The 
spin-rotation coupling term is exponentiated as a chem- 
ical potential [loj . The lattice quark action is 



Sf i){xi) 



«|(l-7^)T,+ + (l + 7^)r, 

(1-7^^)T,+ + (1+7^^)T,_ 
(l-7^)T,+ + (l+7^)r,_ 



(18) 



(1 - 7^) cxp ( mil— ) Tr+ 
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(1 + 7"^) cxp ( — iafi-^ ) Ti-- 



with = Uf^{xi)5xi+ii,x2 and = Ul{x2)5xi-fL,x2- 
The bare hopping parameter is k = l/{2am -\- 8). In 
the continuum limit a — >■ 0, the lattice actions (jlSp and 
([T5)) correspond to the continuum actions ([7]) and ([TU| . 
respectively. 

In the Minkowskian rotation, the angular velocity is 
replaced as $7 — ifi. In the gluon action ([T5|) . the 0(17) 
terms become pure imaginary numbers. In the quark 
action (jlSp . the orbit-rotation coupling term becomes 
like an imaginary hopping term and the spin-rotation 
coupling term becomes like a chemical potential. Since 
both of the gluon and quark actions become complex, 
the Monte Carlo simulation severely suffers from the sign 
problem. On the other hand, in the Euclidean rotation, 
the gluon and quark actions are real, and thus there is 
no sign problem. 

We have formulated the hypercubic lattice, which is 
commonly used in most lattice simulations. It is possible 
to formulate the cylindrical lattice in the cylindrical co- 
ordinate — (r, 0, z, r). An advantage of the cylindrical 
lattice is better rotational symmetry around the rotation 
axis. However, since the action includes a singular metric 
factor l/r, the region around the rotation axis must be 
removed to avoid this apparent singularity. 

The angular velocity $7 might affect the renormaliza- 
tion, e.g., the physical scale. This correction cannot be 
neglected when the angular velocity is large. However, it 
is not trivial how to determine the physical scale in rotat- 
ing frames. Moreover, since the Lorentz symmetry and 
translational invariance are broken, the isotropy and the 
coordinate independence are no longer assured at the full 



quantum level. (This is similar to the anisotropic lattice 
[m and the coordinate-dependent lattice coupling [l^.l 
In the following numerical simulation, we restrict the an- 
gular velocity only to small values, and do not discuss 
the problem of the renormalization. 

Simulation. — We performed the quenched SU(3) 
Monte Carlo simulation. The lattice size is x Ny x 
X A^T- = 13 X 13 X 12 X 12. The range of the x-y plane 
is a; = [—6a, 6a] and y = [—6a, 6a], and the position of 
the rotation axis is (x, y) = (0,0). We set the bare lat- 
tice coupling j3 = 5.9 and the bare hopping parameter is 
K — 0.1583, where the lattice spacing is a ~ 0.10 fm and 



the meson mass ratio is vn-^jmp ~ 0.59 13 1 



We analyze the angular momentum of the rotating 
QCD vacuum. Rotation induces a finite vacuum ex- 
pectation value of the angular momentum operator. To 
understand the reason, let us recall the energy of a ro- 
tating classical particle. In a rotating frame, the en- 
ergy is = mr^^^/2 = mr^{()rcst — ^7)^/2, and it has 
a minimum at ^rcst — The lowest energy state 

in the rotating frame has a finite angular momentum 
L = rnr^^rost = rnr^fl in the rest frame. Similarly, in 
QCD, we can observe the rotating state with a finite an- 
gular momentum by generating the lowest energy state, 
i.e., the vacuum, in a rotating frame. 

We calculate the vacuum expectation value of the an- 
gular momentum density operator. Since we want to 
know the angular momentum observed from the rest 
frame, the operator is defined in the rest frame. 
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(19) 



where £ is the Lagrangian density. We calculate the ex- 
pectation value of this operator in the rotating vacuum. 



(20) 



The angular momentum density operator (|19p coincides 
with the conserved Noether current in the flat space-time. 
The gluon angular momentum density is 



Jg 



1 



IxFyzF^r 



'^yFxzFzr] 



(21) 



The fermion angular momentum density is decomposed 
into the orbital and spin angular momentum densities. 



Jf — Jfl + Jfs, 
Jfl = {-i'Y{xDy-yD,)i^), 



Jfs 



-i^-f^^^ 



(22) 
(23) 

(24) 



We discretize these operators in the same way as the 
lattice actions (|15p and ([TS]) . 
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FIG. 2: Angular momentum density J along the x axis with 
the angular velocity aO, = 0.06. The solid curves are quadratic 
fitting functions. 
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FIG. 3: Angular momentum density J at (a;, y) = (2a, 0) as a 
function of the angular velocity $1. The solid curves are linear 
fitting functions. 



In Fig. [21 we show the angular momentum density 
along the x axis (y = 0). The angular velocity is fixed at a 
nonzero value af2 = 0.06. As for Jq and Jfl, the angular 
momentum density is a quadratic function of the distance 
from the rotation axis. The spin angular momentum den- 
sity Jps is small but nonzero, and it is independent of 
the distance. In Fig. [31 we show the angular momentum 
density measured at a certain point, [x, y) = (2a, 0), as a 
function of the angular velocity 51. The angular momen- 
tum density increases linearly. From fitting the data, 

Jg = (0.94±0.01)a-^ X r^ll, (25) 
Jfl = (0.60±0.01)a-^ X r^ll, (26) 
Jfs = (0.17±0.01)a-2 X fl. (27) 



The coefficient in front of Q. is interpreted as the moment 
of inertia of the constituent in the QCD vacuum. The 
functional form of J a and Jp^ can be intuitively under- 
stood from the angular momentum of a classical particle, 
L = mr^fl. The numerical coefficients of Jq and Jfl are 
interpreted as the inertial mass densities of glueballs and 
quark-antiquark pairs, respectively. The r-independence 
of Jfs is a plausible result since the spin is an intrin- 
sic angular momentum. Note that these coefficients are 
unrcnormalizcd and they depend on the renormalization 
scale and the quark mass. 

Summary. — We have formulated lattice QCD in rotat- 
ing frames. We have carried out its first Monte Carlo 
simulation to analyze the angular momentum of the ro- 
tating QCD vacuum. At least in the case of the Euclidean 
rotation, we can implement this framework without tech- 
nical difficulty. By using this framework, we can study 
rotating hadrons, rotating QCD matter, and many other 
rotational phenomena from first principles. 
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